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Abstract 


The vortical flow over a delta wing contributes an important part of the lift - the so 
called nonlinear vortex lift. Controlling this vortical flow with its favorable influence 
would enhance the maneuverability and control of the aircraft. Several experiments 
have revealed that span wise and tangential leading edge blowing can be applied as 
a means of controlling the position and strength of the leading edge vortices. The 
present research studies these effects by numerical sedation of the three-dimensional 
Navier-Stokes equations. The leading edge jet is simulated by defining a permeable 
boundary, corresponding to the jet slot, where suitable boundary conditions are im- 
plemented. Numerical results are shown to compare favorably with experimental 
measurements. It is found that the use of spanwise leading edge blowing at moderate 
angle of attack magnifies the rise and strength of the leading edge vortices, and moves 
the vortex cores outboard and upward. The increase in lift primarily comes from the 
greater nonlinear vortex lift. However, spanwise Mowing causes earlier vortex break- 
down, thus decreasing the stall angle. The effects of tangential blowing at low to 
moderate angles of attack tend to reduce the pressure peaks associated with leading 
edge vortices and to increase the suction peak around the b— edge, such that the 
integrated value of the surface pressure remains about the same. Tangential leading 
edge Mowing in post-stall conditions is shown to re-establish vortical flow and delay 
vortex bursting, thus increasing C V — and stall angle. 
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Chapter 1 
Introduction 


1.1 Motivation 

An improved understanding of the aerodynamics of slender, highly swept wings over 
a broad range of angle of attack is necessary fix the design of aircraft haying greater 
maneuverability and aerodynamic control. The low aspect ratio delta planfonn wing 
is typical of that used for highly maneuverable aircraft; the flow about such a wing is 
characterized by the presence of large scal e , organized vortices on the leeward side of 
the wing which occur at moderate to high angles of attack. The formation of these 
vortices is responsible for the favorable “nonlinear" effects on the lift and may be 
used to advantage in both maneuvering and configurations. How e v er, at the 

very high angles of attack attained during post-stall maneuvers, vortex breakdown 
can occur; ie., the vortices become unstable and unsteady, resulting in undesirable 
forces and moments acting on the ving. 

The flow pattern about a delta wing can, in general, be c l as si fied into four 
regimes as a function of angle of attack, (table 1.1). 
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(I) At very km angle* of attack, the flow remains attached and vortex-free. The 
lift is linear with angle of attack. 

(II) At moderate to high angles of attack, the boundary layer on the lower surface 
flows outward and separates at the leading edge resulting in a free shear layer 
emanating from the edge. This free shear layer raOs up in a spiral fa s h i on 
to forma pair of vortices on the leeward side of the wing, as shown in figure 
1.1 . These vortices are characterized by being both stable and symmetric. A 
noticeable difference with regime (I), as illustrated in figure 1.2, is that the 
lift increases naolinearly as angle of attack increases, due to the low pr es su r e 
on the wing upper surface induced by the leading edge vortices. This is the 
so-called “nonlinear lift”. 

(III) At very large angles of attack, breakdown of the vortex structure occ u r s near 
the trailing edge and may be accompanied by an asymmetric and unsteady 
flow pattern. A distinct indica t ion of the vortex breakdown is an abrupt 
increase in the cross sectional area and turbulence level of the vortex co re s . 
With further increase in angle of attack, the vortex breakdown progresses 
upstream and eventually the vortex system is totally burst. The adverse 
effects of both vortex breakdown and asymmetric vortices result in large, 
unfavorable changes in forces and moments on the wing. 

(IV) At extremely large angles of attack (up to 90*), the flow on the lee ward side is 
characterized by an unsteady wake and periodic vortex shedding. As a result, 
the forces acting on the wing are very unpredictable. 
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Highly maneuverable aircraft operate primarily in the flow regimes (II) and (III), 
where the “nonlinear” vortex lift can be used to aerodynamic advantage. Smce the 
vortex flow for a delta wing plays such an important role in determining lift in these 
flow regimes, it is desirable to control these vortices with their favorable influences 
and thereby enhance maneuverability and control. 

This leads to the motivation o I present work which is to study leading edge 
blowing, either spanwise or tangential to the wing surface, as a means of controlling 
the influence at the leading edge vortices. There are two aspects of interest : one 
is to enhance the vortical lift at low to moderate angles of attack (in flow regime 
II); the other is to control and stabilize these leading edge vortices at high angles 
of attack to delay vortex breakdown (in flow regime HI). 

Since the edge is the region where the vortices are originated, it is logical 

to have a direct control over the vortical flow at its source For other blowing 
concepts, interested readers are r e ferr e d to spanwise blowing [3] - [10] that utilizes 
a discrete jet in a variety of blowing directions, and direct vortex core blowing [11] 
which is used to control vortex breakdown. 


CHAPTER 1. INTRODUCTION 


4 


1.2 Previous Work 

1.2.1 Experimental Evidence 

Span wiia Leading Edge Blowing 

Preliminary experiments by Trebble [12] on a 70* swept delta wing in the angle- 
of-attack (a) range of 0* to 16.4* have shown that lift gains (increase in Ct at a 
given incidence) axe obtainable by ejecting a high-momentum jet in the spanwise 
direction from the leading edge of the wing, as depicted in figure 1.3 . The lift was 
measured by a force-balance and results demonstrated that there was no appreciable 
Reynolds number effect in the range of 0.35 x 10* to 1.05 x 10*. five different jet 
nozzle shapes were tested and the optimum nozzle shape proved to be a tapered 
slot which was expected to give a momentum distribution ap p r opriate to preserve 
conical flow. 

Alexander [13, 14] conducted a aeries of experiments on a 70* swept delta wing 
and a cropped delta wing with leading edge blowing for a = 0* to 25*. The results 
revealed that blowing increases the size and strength of the leading edge vortices 
and moves the cores outward and upward. The increase in lift is mainly due to an 
increase in the nonlinear contribution. In addition, the suppression of the secondary 
separation due to the entrainment of the leading edge jet was also o b ser v ed in the 
exper im ent 

The same concept was applied to a straight wing and analysed theoretically and 
experimentally by Tavdla et si [15,16, 17] Although the dependence of the lift on 
the vortical flow for a delta wing and a straight wing is quite different, the outward 
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displacement of tip vortices by the jet and the effective increment of the wing span 
are observed in both cases. 


Tangential Leading Edge Blowing 

This blowing concept modifies the vortical flow by controlling the sepa r ation of the 
cross-flow boundary layer with a thin, high- momentum jet ejected tangentially to 
the leading edge surface. Experimental work on tangential leading edge blowing is 
more recent, and these investigations were carried out by Wood ef si [18, 19] These 
tests were conducted on a 60* -sweep conical delta wing with a rounded Wting edge, 
as depicted in figure 1.4 . Experimental results show that the use of tangential 
blowing at post-stall angles of attack (tested up to a =* 60*) re-establishes vortical 
flow and delays vortex breakdown, thus increasing the Cr and the stall angle. 
At pre-stall angles of attack, tangential leading edge blowing tends to reduce the 
suction peaks at the wing surface associated with the vortices, while increasing the 
suction peak around the leading edge such that the integrated value of the surface 
pressure remains about the same. The implication of this pre-stall characteristic of 
the lift response to blowing is that, when blowing is applied at a given post-stall 

angle, there is an optimum blowing coefficient to reach Cr Once this maximum 

Cr is reached, additional lift can not be obtained with further increase in Mowing 

intensity. In other words, as the vortical flow is re-established due to blowing in 
post-stall conditions, the flow topology resembles that at lower angles of attack, but 
the lift is greater. 




Span wise Leading Edge Blowing 

Barsby [20,21] analyzed the flow field over a delta wing with spanwise blowing using 
the assumption of a conical flow field. He adapted the vortex-sheet model, applied 
by Mangier and Smith [22] to the study of leading edge separation, combined with 
the thin- jet model introduced by Spence [23] in the study of the jet flap. The jet is 
modeled as an infinitely thin vortex layer which carries a pressure difference between 
its upper and lower surfaces. This p r es sur e jump across the sheet is proportional 
to the product of the momentum flux of the jet and the local curvature of the sheet 
and is related to the distribution of vortidty on the sheet. 

In Barsby’s solution procedure, the vortex sheet trace in the cross-flow plane is 
divided into vortex segments, and a concentrated vortex core is introduced which is 
joined to the last segment by a cut. Three separate nested iterations are performed : 

(1) The pressure jump across the jet-vortex sheet is calculated at the center of 
each segment using the initial approximation for the shape of the sheet. 

(2) The vortex position is adjusted by satisfying the force-free condition on the 
vortex segments, core and cut. 

(3) The sheet shape is adjusted to satisfy the velocity tangency condition. 

The procedure is repeated until the vortex sheet shape no longer varies. Blowing 
is shown to increase the lift on the wing and to increase the circulation about the 
vortex core, while displacing it upward and outboard. 
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Although Baabjr’a result* closely re s em ble the sctusl flow field, general agree- 
ment b et we en theoretical predictions and experimental results for subsonic flows is 
satisfactory only for the front portion of a delta wing. This limitation is due to the 
absence of the trailing edge in the conical assumption. 

Recently, theoretical ex p ressi ons for the effects of span wise leading edge blowing 
were derived by Tavella [24]. In his analysis, the separated flow on the delta wing 
is r e pr e sen ted by a model consisting of a pair of concentrated vortices c on nected 
to the leading edges by straight vortex sheets, originally proposed by Brown and 
Michael [25]. The blowing problem is solved implicitly in the cross-flow plane by 
requiring the singularity system to be force-free in the presence of the jet momen- 
tum. T>x» predicted quantities d r»mwV»hU agreement with that m tTw» 
experiments [12]. However, the lack of a detailed structure of the flow Add in this 
analysis limits the insight that can be gained in the physics of spanwise blowing. 

Barsby’s and Tavella ’s analyses assume conical flow and therefore cannot account 
for the vortex breakdown phenomenon. 

Tangential Leading Edge Blowing 

Mourtoa [26] studied the control of vortical separation on a circular cone by tangen- 
tial blowing. This model utilises an iterative p roce dur e by matching the solutions 
of the viscous and the inviscid flow fields. The single line vortex model [27] was 
used for the inviscid solution, which is similar to that of Brown and Mi chae l [25]. 
The integral form of the cross-flow boundary layer equation is solved to determine 
the location of the cross-flow se p aration point. This iteration continues until the 
pr es s ure on the upper and lower separation points are matched. 
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The controlling of the separation point through tangential blowing is demon- 
strated in Mourtoa’ analysis as a means of influencing the strength and position of 
the vortices at low to moderate angles of attack. The predicted separation lines 
agree well with experiments [18, 19]. How e v er, due to the simplicity of the invisdd 
model, the detailed flow field agrees only qualitatively with experiments [18, 19] (*.e. 
vortex loca t ion and strength). In particular, the a b sence of trailing edge condition 
precludes the use of this method to study vortex breakdown. 


1.3 Present Approach 

The objective of the present research is to study the effects of spanwise and tangen- 
tial leading edge blowing on the vortical flow field over delta wings. The approach 
is to solve this three-dimensional problem using Computational Fluid Dynamics 
(CFD) techniques. The 1— /ting edge jet is simulated by defining a permeable 
boundary, c orre s po n ding to the jet slot, where suitable boundary conditions are 
implemented. A single-grid is used for the spanwise blowing calculations, whereas 
a multiple-grid zonal type of approach is required in the tangential blowing cases 
to accomodate the jet -slot geometry. A turbulence model for a wall jet developed 
by Roberts [28] is incorpor a ted to accurately predict the location of the separa- 
tion point in the tangential blowing calculations. Aided by available experimental 
evidence, the structure of the flow field is studied and analyzed through detailed 
graphical repres e n tations. The pressure field, including the surfoce p re ssu re distri- 
bution on the wing, is computed and the variation of lift with angle of attack is de- 
termined, with and without blowing. These results are compared with the previo us 
experimental results of Trebble [12] (for spanwise blowing) and Wood et al [18, 19] 
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(for tangential bkming). Such nnm*nf1 experimentation not only enables one to 
gain physical insight in the flow field structure, but also reduces the number of 
configurations or models to be tested during preliminary design. 

This work is divided into five chapters. The g overn i ng equa t ions sad numerical 
algorithms are presented in Chapter 2. Chapter 3 describes the implications of the 
spanwise lead ng edge blowing concept an vortical lift and its physical interpretation. 
Chapter 4 describes the effects of tangential leading edge blowing on the vortical flow 
field and its applications to post-stall condit i ons. Finally, Chapter 5 summarizes 
the results, presents conclusions and dismisses r ecommen da t ions for future work. 
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Regime 

Angles of Attack 

Flow Characteristics 

Delta Wings Forces 


moderate 


IV 

extreme 


(up to 90*) 




symmetric 
vortical flow 


Vortex asymmetry and/or 
vortex breakdown 


[if.! 11-' M 


ear 


Nonlinear 


Table 1.1: Flow regimes as a function of angle of attack. 


Flow Regime 


nonlinear part 


linear part 


Figure 1.2: Schematic repr 
(Hummel, Ref. 2). 


station of the lift characteristics of a delta wing 














Figure 1.3: Schematic view of a delta wing with spanwise leading edge blowing. 





Chapter 2 


Governing Equations and 
Numerical Method 


2.1 Governing Equations and Approximations 

The three-dimensional Navier-Stokes equations are used in this problem to describe 
the physical phenomena of the flow. The equa t ions are solved with the finite- 
difference algorithm developed by Fujii and Obayashi [29,30,31]. To solve the 
Navier-Stokes equations through a finite-difference or a finite-volume approach, the 
equations are generally put into nondimensional, con s er v ation-law form and the 
computation is carried out in a generalized coordinate system, using the thin-layer 
approximations for high Reynolds number flows. The governing equations and the 
numerical algorithm will be discussed briefly. More detailed information is found in 
references [29,30,31]. 
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2.1.1 The Navier-Stokea Equations 

Under the assumptions of no body force and no external heat addition, the Navier- 
Stokea equations in nondimenaional Cartesian coordinates can be written as (see 
Ret [32]) 

0Q_ OE OF OG OE, OF. OG, . 

9t + 0x + dp* Ox Ox + Ov + Ox 

where the vector of conserved quantities Q is 


<?= pv 


the Euler flux vectors E, F, G are defined as 


P“* + P 


«(* + p) 


, F - 


+ P , G = 


»(« + p) 


pw 3 +p 
t»(e + p) 


and the viscous flux vectors E„ F„ G, are given by 


E, * Re” 1 Twm , F, - Re- 1 r „ , G, = RT X \ r. 
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with 


r„ = A (u, + v, + w,) + 2pu, 

T„ as A (U, + t> f + W t ) + 2fiV t 
T„ a A (li, + 0, + U» t ) + 2flW, 

T *t - = /i (u, + v M ) 

r„ sa T„aa/l(tl, + t»,) 

T w . = r n = n(v, + w t ) (2.5) 

p x as 7«Pr _1 3 r «, + UT„ + ur*, + urr„ 

/?, as 7/cPr _1 3 f e t + ur,* + vr„ + urr y , 

P, = 7 KPr~ l d,t % + ur„ + + icr„ 

e, s= ep -1 — 0.5 (u J + v 2 4- w 2 ) 

The density is nondimensionalized by the free stream density the Cartesian 
velocity components u, v and w are normalized by the free stream speed of sound 
Ooq, and the total energy per unit volume e, is scaled by Poooie- Pressure is related 
to the conservative flow variables, Q, through the equation of state for a perfect gas 

P = (7-l)[«-|(« J + w S + ®*)] (2-6) 

where 7 is the ratio of specific heats. A Newtonian fluid is considered, with the 
coefficient of bulk viscosity A obtained from the Stokes’ hypothesis, A = — J/i. The 
dynamic viscosity ft is nondimensionalized by Moo, * is the coefficient of thermal 
conductivity, Re is the Reynolds number and Pr is the Prandtl number. 


A 



m.- 
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2.1.2 Coordinate Transformation 

In the numerical computation, the governing equations are usually transformed to a 
general body-conforming curvilinear coordinate system. Thia males the formulation 
of the code general and independent of body geometry, while increasing numerical 
accuracy (for the approximation that will be introduced in section 2.1.3 ) and fa- 
cilitating the handling of boundary conditions. Tins is achieved by the following 
transformation : 


t = £(*, y,*,t) 

n = q(x,y,z,t) (2-7) 

C = 

Thus, the governing equations are transformed from the physical domain (x, y, z) to 
the computational domain (f,q,C)- Expanding Eq. (2.7) via the chain rule yields, 

§t 1 ft Vt Ci 37 

* = 0 U * 0 * ( 2 . 8 ) 

£ • 0 i» (i & 

. & . , 0 & V* Ci . . ft . 

or conversely, 


1 *r »r it 

0 *< y t z e 
0 x, y, z, 
0 x c y c z ( 




» 
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The transformation Jacobian is given by 

Using the transformation in Eq. (2.9) and dividing by the transformation Jaco- 
bian, the Navier-Stokes equations can be written in strong conservation-law form 
[33] in a generalised coordinate system ((,q, Q as follows 


qq d(P-P.) d(6-&.) 

&F + at + zh + ac 


( 2 . 10 ) 


where 


0* J" 1 


/m , 


E^J- 1 


puU + tsp 
pvU + i t p 
pwU + i.P 

U(e + p)-£,p 


( 2 . 11 ) 


P=J~ l 


pV 

puV + ij,p 
pvV + TIrf 
pwV + q.p 

v(*+p)-v*p 


C = J~ X 


pW 

puW + CsP 
pvW + C,P 
pwW + C *p 
H^Ce+pJ-CiP 


( 2 . 12 ) 


The contravariant velocity components U, V, W are defined as 


U = 6 + &« + & + 

V -Tft + 1 /,« + 1J 9 V + Tf t W 

W = o + c*« + c*» + c® 


(2.13) 
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The viscous flux rectors ere given by 


(i r M + + 6 r a 

= J-'Rz-' e.r^ + « t r„ + 6v 

(i r a + + (< r n 

CsA + + GA 


(2.14) 


9t r n + Vq + 7»*V» 

F. = J~ l Re~ l i) a T w + 7,t w + 7,r„ 
7« T *» + r 7* r qf + T 7* T »» 
»7» A + 7»A + »7* Pa 


(2.15) 


C.’w + C,r 0 + C,T„ 

G. - <,r„ + C» r w + G«W 

G^- + G r *» + G r « 
GA + GA + GA 


(2.16) 


In the viscous terms, the Cartesian derivatives are transformed to the generalised 
coordinatee using the chain rule, Le. 


tt. = G«< + »?««s + G«< 


Use of the inversion relationship between Eqs. (2.8) and (2.9) leads to the metric 





Tr' 


4 


3 
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identities as follows 


=■ J(x c z,-*,z c ), 

6 »^(*syc -*<¥<»). 
rh = J(v<*t- y<*c), 
= -*<*<). 
»?*“^(*<y<-*<y<). 


C. - ^ (y<*, - y^t) 

-*<*,) 

C. * J (*<», - *,»<) 

6 = -*r4*-yr&-*r 4* 
•?» = ”*»•)* Vf* 7 , ~ *r* 7 » 

C« = -*rCr - yrCf - *rC» 


(2.18) 


The inverse of the transformation Jacobian is 


J- 1 - det 


*i *, *< 

y< Vs v< 

L *t *< *c 


= *«y^ z c + *<y<*» + **vc*< - *«yc z s - *,y<*c - x <yv z t (2-i9) 





r 


2.1.3 Thin-Layer Approximation 

For high Reynolds number flows, the viscous effects aze confined to thin vortidty 
layers on the walls and in the field. Resolution of the normal component of velocity 
gradients, and hence of boundary layers, demand high grid point density in the 
normal (£) direction near wall surfaces. Due to memory and time constraints this 
makes it impossible to have enough grid density in the streamwise (£) and circum- 
ferential (if) directions to resolve the corresponding derivative components in the 
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viscous terms. From this point ai view it is appropriate to use the thin-layer form 
of the Navier-Stokes equations. 

In the thin layer approach only the (-component (normal to the wall) derivatives 
in the viscous terms are retained. The thin-layer approximation can be summarised 
as follows. 


dt . 

ae 


0 


dK 

dn 


o 


and the Cartesian derivatives for viscous terms are approximated by 


( 2 . 20 ) 


h 


6 C » 


0 

h 

= 

(* Vw (f 


0 

. h . 


, (» *?• C* , 


. k . 


where ^ and Jjj have been set to zero. 

The full Navier-Stokes equations (2.10) then simplify to 


(2.21) 


dQ dt dt . d& B .dt 

dr + dt + dn + dc ac 


( 2 . 22 ) 


where 


0 


$ = J- 1 


+ (M/3)mjO 
nm x v c + (m/ 3) mj(, 
+ (M/3)m,C, 


ftmim 3 + (ft/ 3 ) m } (ij m v + t} t v + 17,1s) 


(2.23) 
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with 


m i =* C» + C? + 

m J = + C#*>< + C* w c (2-24) 

m, - \ (u J + * J + w J ) < + nPr~ l (7 - 1 )’ 1 (a*) ( 

Eq. (2.22) represents the thin-layer Navier-Stokes equations in a generalized coor- 
dinate system. 


2.2 Numerical Algorithm 

The numerical method used in this computation is the LU-ADI factorization algo- 
rithm proposed by Fujii and Obayashi [29]. This is an implicit, noniterative scheme 
which is first-order accurate in time and fourth-order accurate in space. The ba- 
sic structure of this method is baaed upon the Beam-Warming algorithm [34] and 
incorporates the concepts of a flux vector splitting technique [35j and diagonally 
dominant factorization[36, 37] as described in the following sections. 

2.2.1 Beam and Warming Algorithm 

The Beam and Wanning [34] (or Briley and McDonald [38]) implicit approximate 
factorization algorithm is commonly used in Computational Fluid Dynamics. Ba- 
sically, this scheme can be either first- or second-order accurate in time, second- 
er fourth-order in space and is noniterative. At each time step, the equations 
are factored spatially which reduces the solution process to three one-dimensional 
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problems. The spatial derivatives are approximated by t h re e -po i nt, second-order 
central differences, which produces block tridiagonal systems for each spatial direc- 
tion. Although the linear stability analysis in r e f e r enc e [34] shows the schemes to 
be unconditionally stable in two-dimensions, stability restrictions are encountered 
in practice due to the nonlinearity of the equations. H o wev er , the Beam- Wanning 
algorithm is unconditionally unstable in three dimensions, the numerical stability 
is achieved by adding artificial dissipation terms to the equations. 

The Beam- Wanning factorization applied to (2.10) can be written as 


(/ + h£ ( A' - Al<) (/ + AM* - Al„) X 

(/ + + d,I c ) aq- = hr 

= -a (Stir + + s ( dr - Re' t s ( s n ) 

~ (D» |{ + -0« U + D, 1c ) 0* (2.25) 

where and superscript n denotes the n 4 * time step. This method is 

in the so-called “delta form”, which is widely used for solving steady-state problems 
since the final solutions do not vary with the left-hand-side operators once steady 
state is reached (*.e. A Q u =* o). In Eq. (2.25), I is the identity matrix and h 
is the time step, 6 is a central finite-difference operator, A and V are forward 
and backward difference operators, respectively. As an example illustrating these 
operators, consider the (-direction 

6 t Q n = [$*(( + A(,q,C)-0*«-A(,q,O]/2A( 

A tQ n = [<5"(( + A(,q,C)-0*(f,Tf,C)] /** 


(2.26) 
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ViQ n - /A4 

The flux Jacobian matrices A", B n and 0 n are obtained from a local linearization 
by taking a Taylor series expansion about 


£r+ l = ir+A n (Q n + l -Q n ) + o(h a ) 

F"* 1 = F* + F , (Q" +, -$ - )+0(h a ) (2.27) 

G”+ l = G n + C n (Q* +t -Q n )+o(h 7 ) 


where the flux Jacobian matrices are given by 


A = 


dE 

dQ 




(2.28) 


Similarly, the viscous coefficient matrix M is obtained from the linearization of the 

A A A A 

viscous vector 5. Expressions for the flux Jacobian matrices. A, B, C as well as the 
viscous coefficient matrix M are given in Appendix A. 


The D, and D, terms are implicit and explicit artificial dissipation terms which 
should be added to the left-hand-side and the right-hand-side, respectively, to main- 
tain numerical stability. The expressions and origins of these numerical dissipation 
terms will be described later. 


2.2.2 The LU-ADI Algorithm 

The procedure for solving the block tridiagonal systems associated with the Beam- 
Warming algorithm [34] requires a large computational time. In order to make the 
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solution scheme mote cost effective, the left hand side of Beam* Wanning scheme 
is mrwliftHrf by introducing the LU-ADI algorithm described in reference(29]. This 
algorithm incorporates a flux-vector splitting technique [35] and a diagonally domi- 
nant factorization [30, 37] such that the Alternate Direction Implicit (ADI) operator 
is decomposed into a product of lower and upper Indiagonal matrices. In this way, 
only scalar inversions are required. 

The flux Jacobian matrices A, & and £ given in Eq. (2.28) have real eigenvalues 
and a complete set of eigenvectors and can then be diagonalized through a similarity 
transformation, 


Ac-i^ATc, A, = i7 l dr,, A< = 17 1 £ r c (2.29) 


The diagonal matrices A t , A„ A< contain the eigenvalues in each spatial directum 


A< = Diag [U, U, U, V + ay/Q + $ + Q, U-a^+^Q] 

A, = Diag [v, V, V, V + + V - + + (2.30) 

A< = Diag [W, W, W, W + a,/<2 + $ + (*, W - a,/<2 + <* + <3 ] 

where U, V and W are the contravariant velocity components and a is the speed of 
sound. The similarity transformation matrices T<, T„ T( and their inverse matrices 
TJT 1 , T^ x , are given in appendix B. 

For the (-direction, the Beam-Warming ADI operator (£) can be rewritten in 
diagonal form [36, 37] as follows 
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t ( m I + kS(A + 

- T(T[ l + hS t [T € A t T{ 1 ) + TtJ-'^JTf 1 
« T<[l + hS t A« + J- l €i6\j] 27 1 (2.31) 

where an approximation is introduced by factoring the eigenvector matrix T t outside 
of the spatial derivative operator <$;. Since the eigenvector matrix T { is a function 
of this modification reduces the time accuracy to at most first-order in time [37]. 
It should be pointed out that the explicit side of the diagonal algorithm which con- 
tains the steady-state finite difference equations remains the same. If the diagonal 
algorithm con ver g es, the steady-state solution will be identical to that obtained 
from the original algorithm. Thus, this modification does not downgrade the final 
steady state solution. 

The central differencing can be decomposed into two-sided differencing by using 
the flux-vector splitting technique [35]. The right hand side of Eq. (2.C1) becomes, 

C t = T { [/+ V t A+ + A t A<] I7‘ (2.32) 


with 




(2.33) 


where contains the positive eigenvalues and contains the negative eigenvalues, 
XT 1 is taken to be the transformation Jacobian at the central point, and e, is the 



t 
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magnitude of the implicit artificial dissipation term, to be discussed in the next 
section. 

Using first-order upwind difiereneing, Eq. (2.32) can be written as 

CfT t [L ( + D f + U t ] (2-34) 

where 

L t = -A^ 7 _ 1 

D i = /+( AJ,-A7,) (2.35) 

U t = 

where the subscript j, is the grid paint index in the £ -direction. Eq. (2.34) can be 
approximated as follows 

L t + D t + U ( = (L t + D t ) DJ 1 (D t + U ( ) + 0 (h 2 ) (2-36) 

where higher order terms are dropped since D t = 0(1) and L(, U t = 0(h). Sub- 
stituting Eq. (2.36) into Eq. (234), the LU factorization for the ADI operator in 
the ( direction becomes 

I + h6 ( A + J~ x eS[J = T ( (L ( + D ( ) D? (D t + U t ) Tf l (237) 

The hlock tridiagonal system has been decomposed into the product of the lower 
and upper scalar bidiagonal systems [L^ + D(J and [D^ 1 (D( + Z7 f )] . 
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ii 
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In order to take the thin-layer viscous tens* into account while using the diagonal 
dominant factorization, the split flux Jacobian matrices 6* are modified aa 

£* « T< (A* ± vl) 77 1 (2.38) 

where 

*'=2ny/ci+G+(:nR*p±Q 

By using a similar procedure for the other operators, this method can be sum- 
marized as follows : 


T< (L t + D e ) D? (D< + U ( ) (L, + £>„) D~ l (A, + tf n ) 

(17*T C ) (L< + D<) D? (D ( + U ( ) T£ l &Q n * RT (2.39) 

The expr e ssi ons for (T^T'T,) , (T^ l T c ) used in Eq. (2.39) and their inverse matrices 
(T^ l T() , (T^ l T n ) can be written in a simpler form, given in Appendix B. In the 
solution process, an inversion in one direction consists of one scalar forward sweep 
and one scalar backward sweep. 


2.2.3 Numerical Dissipation Model 

The fourth-order dissipation model has been widely used with central differences 
in Navier-Stokes computations. This constant dissipation model produces an unde- 
sirable oscillatory solution near discontinuities such as shock waves. On the other 
hand, second -order dissipation produced by proper upwinding works much better 


s 




1 



JLt- 
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than the fourth-order dissipation model in regions near discontinuities. Such dis- 
sipation, however, is not desirable in the rest cl the flow field since it reduces the 
spatial accuracy to first-ocder. 


The numerical dissipation model used here is a non-linear, second- and fourth- 
order mixed type whose magnitude is controlled by a local flux limiter functi o n 
baaed on a simplified version of the Total Variation and Diminishing (TVD) [39] 
method. In doing so, numerical stability is maintained while the flow physics is 
captured without an excessive amount of artificial dissipation, as compared to the 
global constant numerical dissipation model. 

For the f - direction, the explicit numerical dissipation terms can be expressed 


D. k <r - V« {* *$,+(*- D V<A< [(i)^ M,] } (*•«) 

where 

d = flux limiter function 

e, * K t hou 

K, = an input constant 
ffk =* spectral radius of A* 

= |ui* + vk, + w*,| + o^i* + kf + k*;k = (, tjctC 

The flux limiter function d has values ranging from 0 to I depending on the smooth- 
ness of the solution. For a relatively smooth solution, for which the signs of forward 
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A 






and backward differencing at a given point of the solution are the same, and their 
magnitudes are similar, d takes the value of zero so that only the fourth-order dis- 
sipation terms are used. On the other hand, when the flux limiter function detects 
an oscillation or a large gradient, d will be set dose to 1 so that the second-order 
dissipation terms are dominant. 

In order to reduce the inflii«w» of the numerical dissipation in the direction 
normal to wall inside of boundary layer, the magnitude of numerical dissipation is 
scaled by the ratio of the local contravariant velocity component to the free stream 
velocity magnitude 


««l< 



(2.41) 


where, for example, V( * (£*, £,). The details of the current flux limiter function 

are given in reference [40]. 


In the implicit operators, it may not be necessary to indude artificial dissipation 
terms since the LU-ADI implicit operator is based on upwinding. Currently, implidt 
numerical dissipation terms are included as shown in Eq. (2.33), whose magnitudes 
are obtained in the same way as that for the one in Eq. (2.40). 


2.2.4 Turbulence Model 

The Baldwin and Lomax Model 


i 

i 

i 

r 

! 


1 

i 
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The concept of an effective viscosity coeffident is used to modd the turbulent mixing 


.ii. ! 
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by introducing the so-called eddy viscosity, fit. The eddy viscosity ie obtained from 
the two-layer algebraic model proposed by Baldwin and Lomax [41]. 

In this model, the eddy viscosity is evaluated in inner and outer layers, 


(dinner » C ^ C 

Mt- 

. (Mi)outer . C > C 


(crossover 


(2.42) 


where { is the normal distance from the wall and frrossoTrr “ the smallest value of 
y at which the inner sad outer values cl fi t are equal 

The eddy viscosity coefficient in the inner layer is b ased on the Prandtl-Van 
Driest expre ss ion (Le. reference [42]), 


(Mt)inner = M 


(2.43) 


where the length scale l is obtained from 


where k = 0.4, C + 
vortidty vector 


I- «[i -«*(-£)] 

> * 26 , and |u>| is the magnitude of the local 


M = |V x V| =* >/(u, - v,) 2 + (o, - w,) J + (w, ~ «.) J 


The expre ss ion for the outer layer is given by 


(Mt)outer — K C* Fwm a* Fta w(C) 
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with K s 0.0168, C«, = 1.6 and where F^-l. is determined as follows 


b« = the smaller of 


Com F m 


CU C— Uh / 


The values of (»«, and are obtained from the function 


F(C) = CM 1-exp ( - ^)] 


where C«* = 0.25 and Ui,f is the difference between the marimnm and the minimum 
total velocity magnitudes in the profile along the (-direction. The minimum total 
velocity magnitude is zero at the walls and is a finite in wakes. The quantity F— , T 
is the maximum value of F(() in the profile along the (-direction, and (-,«, is the 
value of ( at which it occurs. The function Fkm(() is the Klebanoff intermittency 
factor given by 


fW()= 1 + 5-5 


(WT 


with C/(itk — 0.3 . 


Once the value of the eddy viscosity, p t , is determined, the total viscosity is 
defined as 




(2.45) 


where pi is the molecular or laminar viscosity. The coefficient of thermal conduc- 
tivity can be obtained as 
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K 


*» + «! = 


Crfil C f Ht 

Pr Pr t 


(2.46) 


The quantity Pr t is the turbulent Praadtl number which is set to 0.9, and the 
lummar Prandtl number Pr is 0.72. 


Turbulence Model for a Wall Jet 


Based on the theory of Roberts [28] for turbulent curved wall jets, an algebraic 
eddy viscosity model, together with the effects of wall curvature, are used in the 
tangential leading edge blowing cases. This model uses self-similar profiles for the 
mean velocity and momentum balance to get the expression for eddy viscosity 


K . 

*t = Tv hpu ' 



(2.47) 


where K = 0.075, k = 0.8814, u m is the maximum velocity, b is the distance to 
the point where u = |u m and ( is the distance normal to wall surface. From the 
existing experimental results (i.e. [43, 44]) it is found 


* = 7 Cm (2.48) 

where the subscript m denotes the location of the velocity. Evaluating 

the constants and making use of Eq. (2.47), Eq. (2.46) can be written as 


Ht = 0.0033 Cm p u m 



(2.49) 


where 
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_cf £ S <<<« 1 

c - 1 ; C — Cm I 

To include the effects of wall curvature, the eddy viscosity is modified as follows 





(2.50) 


where V e is the cross-flow velocity magnitude and R is the radius of curvature of 
the walL This model is applied from the jet exit up to the point where the wall jet 
separates from the surface. 



•* 


1 





Chapter 3 

Spanwise Blowing 


3.1 Introduction 

As discussed in Chapter 1, the vortical flow field asso ci a t ed with flow separation 
from the edge of a delta wing plays an important role in determining the lift. 

It is desirable to control this vortex flow and if possible to enhance its contribution. 
A promising method of control is to eject a thin, high-momentum jet in the spanwise 
direction along the leading edge of the delta wing. This jet sheet displaces the flow 
outboard and rolls up in a similar manner to the rolling-up of the free vortex sheet. 
As a result, higher suction peaks on the wing surface are induced by the str onger 
leading edge vortices which produce more lift. 

Since spanwise blowing increases the strength of the leading edge vortices at 
a given incidence (as will be shown in the later sec t ions,) the advene pressure 
gradient along the vortex core becomes more pronounced near the trailing edge. 
Thus, blowing reduces the angle of attack at which vortex breakdown occurs. For 
this reason, spanwise blowing would be primarily of interest as a lift device, or a 

35 
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roll control device, at low to moderate angles ai attack. Consequently, emphasis 
in this chapter will focus on this angle-cf-attack range. The advene effects on the 
vortex flow due to spanwise blowing in the p resen ce at vortex breakdown will be 
illustrated briefly in section 3.4.2. 

The leading edge Jet is simulated by defining a permeable boundary, corre- 
sponding to the jet slot, where suitable boundary cond iti ons axe implemented. The 
simulations start with the no-blowing case, which will be used to illustrate the 
main features at delta wing aerodynamics and as a base line for comparisons with 
blowing cases. The blowing conditions will be chosen from cases for which experi- 
mental results [12] are available, and will also serve to validate the jet-exit boundary 
conditions. 


3.2 Grid Generation 

The three-dimensional mesh about a delta wing is generated through three steps. 
The first step is to generate the surface-grid that describes the geometry of the 
wing. This is done by distributing the grid points on a two-dimensional section at 
each cross-flow plane, once the ehocdwise stations axe identified. An exponential 
stretching function is used to cluster the grid paints near the edge, where 

the jet slot is located. 

The next step is to generate a three-dimensional mesh nring a two-dimensional 
hyperbolic grid-generation method [45,46] for each ehocdwise station. This method 
uses the surface-grid points as a set of initial data and generates the grid by marching 
away from the body. An inversion of a 2-by-2 block tridiagonal system is required 
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once there are two unknown (p,s) at each grid point under the eontrain of 
orthogonality and poritire cell volume. As a result of using this hyperbolic grid* 
generation method, the outer boundary in the stxeamwise d ir ection is generally not 
uniform, as shown in figure 3.1. This is because the outer boundary is determined 
according to an initial distribution of points on the body, while the local spaa of a 
delta wing in creas e s toward the trailing edge. 

The final step is to redistribute the grid points along each of the radial curves 
by using a cubic spline such that the outer boundary d imen s i ons can be controlled 
and the orthogonality condition dose to the body be satisfied An expon e n ti al 
stretching function is used to duster the grid paints along each radial curve near 
the wall, appropriate to the thin layer approximation Figure 3.2 shows a dose-up 
view of a typical cross section of a delta wing, where bilateral symmetry is in osed 
to reduce the computational effort. 

The final grid arrangement, as shown in figure 3-3, is cylindrical and corresponds 
to an H-type topology in the chord wise direction and an O-type topology in the 
span wise direction. The grid extends 1 root chord upstream, 1 chord downstream, 2 
chords above and laterally, and 1.5 chord below, and results in a total number of 55x 
65 x 33 » 118,000 grid points. The number of grid paints in each direction is chosen 
judiciously to capture the jet effects, while maintaining the memory requirements 
within the 4 million 64-bit words of the Cray X-MP at Ames Research Center. 
Regarding the influence of the over-all dimen s ion s of the c omput a tion al domain, 
numerical experimentation shows that increasing the dimensions in all directions 
by one half of a chord length, while keeping the number of grid points constant, 
changes the lift coefficient by 0.5% (for AR » 1.4, a = 15*). 
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The global re f e r en ce frame is e xp re ss ed in terms of a body-fitted cunrilinear 
system ({, 17 , () where { is in the chord wise direction, 17 is in the circumferential 
direction, and (is in the radial direction such that a right-handed coordinate system 
is formed. 


3.3 Boundary and Jet-Exit Conditions 

Far a steady viscous flow, no-slip boundary conditions are imposed at the wall, 
where U * V = W — 0. The pr e s sure on the body surface is obtained from the 
normal momentum equation. 


or 


[a t i;+a^ , +$6] r 



(3.1) 


(U. + U, + Us) Pi + (»?«C* + TsC, + IsCs) Pn + (C + C ? + C s) P< 

■ -pu (C.«< + + (,»<) - (C.U, + (,*„ + C.w,) (3-2) 

- Pny/G + G + G 

* 0 

where P, is the component of the pressure gradient normal to the body surface. 
Eq. (3.2) is solved at the surface using central second -order accurate differences 
in ( and 17 , and one-sided second-order accurate differences in (. Once the surface 
pressure distribution is determined, the total energy per unit volume, e, is calculated 
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according to Eq. (2.6). Adiabatic or constant temperature walk are used such that 
the density value on the wall k extrapolated from the nearest paint normal to the 
walL These boundary conditions are updated explicitly at each time step. 

At the upstream and circumferential outer boundaries, free stream values are 
specified and remain fixed. Extrapolated values are used at the downstream bound- 
ary except for e, which k calculated based on the free stream static pre ssur e for the 
present subsonic calculations. Symmetry conditions are imposed on the symmetry 
plane. 


Regarding the boundary conditions at the jet exit, it k necessary to specify two 
thermodynamic variables and one kinematic variable at the jet exit plane in order 
to obtain the Q vector. The blowing momentum coefficient ( C„ ) is introduced as 
a nondimens ion&l parameter that indicates the intensity of blowing. 




I f A, P»tV»tV*f*.dA } 
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(3.3) 


where the numerator k the total momentum of the jet, e, k the unit vector in 
the spanwise direction, A } k the projected area of the jet slot perpendicular to the 
span wise direction, k the free stream dynamic pressure and k the wing 

reference area. The jet-exit boundary co nd itions are summarized as follows: 

(1) At the jet boundaries the local static pr e ssu re k used, and the pre ssu re across 
the thin jet k assumed to vary linearly. 

(2) The free stream value of density k specified at the jet exit. 

(3) The jet-exit velocity distribution k prescribed for convenience as an elliptical 
profile whose maximum value k determined from a specified C 0 . 
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As mentioned in the introduction (section 1.2.1), it is desirable to maintain con- 
ical blowing conditions over most of the wing. This requires a linearly varying jet 
momentum flux. This conical blowing condition is approximated through a linearly 
varying jet-slot width and a constant maximum velocity. This distribution of blow- 
ing momentum would be exactly linear if the jet exit velocity profile were uniform. 
The elliptical distribution was chosen, h o we v e r , for better numerical conditioning. 
The analytical calculation of the jet momentum is given in appendix C. 


3.4 Results and Physical Interpretation 

3.4.1 No Blowing 

To understand the effects of leading edge blowing an the vortical flow of a delta 
wing, a no-blowing case at a as 10* is first simulated and will be used for compa- 
rison with the blowing cases and with experimental results [12]. A 70*-sweep delta 
wing is used in the computation with the free stream conditions of Af*, = 0.3 and 
Re*. — 1.3 x 10®. Between the apex and 95% of the root chord ( cq ), the 
delta wing geometry is generated numerically with conical symmetry; that is, with 
linearly increasing thickness of up to 3.5% of cq. The remaining 5 % of the root 
chord is closed by an elliptical afterbody. The leading edge shape is modeled as an 
ellipse of 5 to 1 eccentricity ratio to simulate a leading edge of small radius. 

Figure 3.4 illustrates the leading edge vortex sheet in the half-space from the 
symmetry plane corresponding to 10 degrees angle of attack . The boundary layer 
on the lower surface separates at the leading edge resulting in a free shear layer 
emanating from the edge. This free shear layer rolls up in a spiral fashion to form 
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the lending edge vortex. The particle trace* shown in the figure were obtained bjr 
“releasing particles’* near the wing surface (velocity is aero an the surface) dose 
to the leading edge and tracing their trajectories through numerical integration o f 
the velocity field. This was done by taring the PLOT3D [47] interactive plotting 
software developed at NASA, Ames Research Center. 

A conical transformation [48] is p e rfo r m ed to better illustrate the size and loca- 
tion of the leading edge vortex in the cross-flow plane. The so-obtained description 
of a conical particle trace represents the rescaling of space curves describing particle 
paths of a three-dimensional conical flow [49]. Such a conical transformation can 
be expressed as follows 


X,Y,Z w* 

Conical variables 3-D variables 

with 

X = * 

Y = l 

x 

z = i 

X 

the value of X is chosen to be 1 for simplicity, since the transformation cat 
to drop out. The unit vectors in the conical projection plane are 


+ 17 + 1? 
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and the conical velocity components are calculated according to 


VVonicml - 

(3.6) 

"'cooicd - ? 

(3.7) 


Figure 3.5 illustrates the vortex traces in the transformed conical projection plane 
at 50 % of the root chord ( x/eo =* 0.5 ) In this representation, size and location of 
the leading edge vortex can be easily identified, as opposed to the three-dimensional 
view (figure 3.4). 

This leading edge vortex induces a suction peak at the wing surface, and con- 
tributes appreciably to the lift. To illustrate this characteristic of a delta wing, 
Figure 3.6 shows the pressure coefficient ( C p ) distribution on the wing surface at 
x/cq as 0.5 for a a 10*. A distinct vortex suction peak is clearly evident on the 
upper surface The adverse pressure gradient (slope of the surface pressure curve) 
outboard of the primary vortex core causes a secondary separation, not shown in 
figure 3.5, creating a local suction peak near the leading edge. The narrow suction 
peak located very near the leading edge results from flow acceleration around the 
leading edge from the lower surface of the wing. This suction peak terminates in 
the vicinity of the primary separation line. 

These flow properties are also inferred from the computed near-surface stream- 
lines. For the same flow conditions, figure 3.7 illustrates a typical example of the 
streamline pattern on the upper surface. Aided by the schematic view of the vortical 
flow pattern [50] over a delta wing, as shown in figure 3.8, the following observa- 
tions can be made about figure 3.7 : The flow drawn into the primary vortex core 
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moves sideways until the secondary sep ar a tion Hue is reached. Beyond this line, 
the flow separates due to viscous effects, creating a secondary separation region. 
The flow reattaches at the secondary attachment line and moves further outward. 
This secondary attachment region ends at the primary separation line, where the 
flow from the lower surface separates and tolls up to form the leading edge vortex. 
These attachment and se p aration lines can be easily identified [51] as follows : the 
attachment line is distinguished by diverging streamlines and the separation line is 
characterized by converging streamlines. 


3.4.2 Span wise Blowing 

The same delta wing geometry is used for the spanwise blowing calculations. The 
jet slot geometry is described as a tapered slot that extends approximately from 3% 
to 94% of the root chord, with a maximum jet width of 5.5% of the local thickness 
of the delta wing. Such a slot geometry is chosen to simulate conical blowing and 
to partially match the model used in the experiment [12]. Seven grid points at each 
streamwise wing-station are used to simulate the leading edge jet. The relative 
location of the jet on the delta wing and the blowing direction are shown in figure 
1.3. 

Figure 3.9 shows the particle traces emanating near the leading edge for C 0 * 0 
and .047 ata> 10*, where the particles are released from the same locations under 
the same free stream conditions. Two observations can be made by comparing 
the particle traces far the no-blowing and blowing cases : blowing results in an 
enlargement as well as an outboard displacement of the rolled- up shear layer. 
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case. This ph en om en on is also revealed by the higher local pressure peak near the 
1 fading edge (C„ * 0.0705), as shown in figure 3.11. 

A comparison of the computed results with Trebble’s [12] experimental data 
for a 70* -sweep delta wing at or a* 10* is shown in figure 3.13. In this figure, the 
increment in Ci due to Mowing at a constant angle of attack is defined as 



Agreement between the predicted and experimental results is good. 

Figure 3.14 illustrates the variation of normal force coefficient with angle of 
attack for the no-blowing and blowing (C* * 0.05) cases. It clearly shows that 
higher force is obtained by blowing in the region of low to moderate angles of attack. 
It also shows that the maximum normal force coefficient for the blowing case occurs 
at a lower angle of attack compared to that in the no-blowing case, indicating that 
blowing reduces the angle of attack at which vortex breakdown occurs. The earlier 
vortex breakdown may be attributed to two factors : 

(i) Since span wise blowing increases the strength of the leading edge vortices at a 
given incidence, the adverse pressure gradient along the vortex core becomes 
more pronounced near the trailing edge. 

(ii) The angle of attack at which vortex breakdown occurs decreases as the as- 
pect ratio increases, as illustrated in refe r ence [52]. Since the jet effectively 
increases the wing span, it causes vortex breakdown at lower angles of attack. 

To demonstrate the advene effect of spanwise blowing in post-stall conditions. 
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£ Figures 3.15 sad 3.16 show the rerer w flow velocity contours for the no-blowing 

< 1 and blowing cases at x/co * .77 far o * 40", where vortex breakdown has o ccurr ed. 

The outer contour is the interface boundary between the forward and reverse flow 
velocity, where the velocity magnitude is zero. The r everse flow bubble is clearly 
enlarged and the maximum negative-u velocity is also greater in the blowing case. 
This indicates that blowing enhances vortex breakdown in the post stall conditions. 
Therefore, span wise leading edge blowing would be primarily of interest as a lift 
device, or a roll control device, at low to moderate angles of attack. 

Since tangential blowing (to be discussed in Chapter 4) is primarily intended 
far high angles of attack, the discussion of the structure of the leading edge vortices 
and vortex breakdown at high angles of attack will be delayed until Chapter 4. 
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Figure 3.1: Outer boundary (side view) generated by a two-dimensional hyperbolic 
grid generation method about a delta wing. 
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Figure 3.4: Particle traces 
wing at a = 10* (C M = 0). 
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Figure 3.5: Particle traces emanating near the leading edge on the conical projection 
plane at x/cq * 0.5, a ■ 10* (C 0 * 0). 



Figure 3.6: Computed spanwise surface pressure coefficient (C,) distribution at 
x/cq = 0.5, o — 10* (C M = 0). 
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(b)C m = .0235 


Figure 3. 10: Effect of span wise blowing on the vortex traces on the conical projection 
plane at x/cg = 0.5, a = 10*, for = 0., .0235, .0470, .094. 
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Figure 3.11: Effect of span wise blowing on surface 

x/c 0 = 0.5, a = 10°, for C„ = 0., .0094, .0235, .0470. 
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■*" Predicted values 

Re * = 1.3 x 10®, a = 10° 

• Experimental data; C L \ NoBUmint = .32 


Blowing coefficient, 


Figure 3.13: Comparison of computed and experimental results for a delta wing of 
AR = 1.4 with leading edge blowing at a = 10*. 
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Angle of Attack (a) 
in degrees 




0 

0.00 

0.00 

10 

0.38 

0.40 

20 

0.88 

0.99 ““ 

30 

1.33 

1.47 

35 

1.50 

1.45 

40 

1.18 

1.30 


Table 3.1: Effects of spanwiee leading blowing on total wing normal force (C„ at 0 
and 0.05). 
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Figure 3.14: Variation of computed normal force coefficient with angle of attack for 
C„a0 and 0.05 . 
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«<M< =0.0 

Utni* = -0.14 
Au = -0.01 


Figure 3.15: Negative-u velocity contour for the no-blowing case at */co = 0.77, 
a = 40*. 






Figure 3.16: Negative-u velocity contour for the blowing case (C 0 — 0.05) at 
x/co = 0.77, a = 40*. 
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Chapter 4 


Tangential Blowing 


4.1 Introduction 

The undesirable forces and moments induced by vortex breakdown and a symm etric 
vortex shedding during post-stall maneuvering adversely affect the controllability 
of the aircraft. Experiments [18, 19] have demon str ated that the vortical flow over 
a rounded leading edge delta wing can be controlled by tangential leading edge 
blowing. This is accomplished by controlling the cross-flow primary separation point 
with a thin, high- m omentum jet issuing tangentially to the leading edge sur&ce. 
As a result, the changes in the location of the cross-flow separation point modify 
the strength, stability and position of the resultant vortex as the flow field adjusts 
itself to an equilibrium state This technique is somewhat analogous to the concept 
of circulation control for rounded trailing edge airfoils, where a thin tangential jet 
near the trailing edge separation point controls the location of the flow separat i on 
point and hence the circulation of the airfoil [53]. 


In this chapter, a rounded leading edge delta wing similar to that used in the 
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experiment! [18, 19] is used to simulate the flow field at pre-stall and post-stall con- 
ditions, with and without tangential blowing. A multiple-grid, sonal approach is 
incorporated to accommodate the jet-slot geometry. Two types of zonal-interface, 
overlap and common boundary, are encountered and the interface boundary condi- 
tions are treated accordingly, as to be described in section 43. The leading edge 
jet is implemented through a set of "permeable” boundary conditions on the body 
surface. For the no blowing cases, the solid-wall boundary conditions are imposed at 
the jet-exit. The algebraic turbulence model of Baldwin and Lomax [41] is utilized 
near the wing surface. The algebraic turbulence model for a wall jet developed by 
Roberts [28], as described in section 2.5, is used in the jet region up to the point 
where the wall jet separates from the wing surface. 

The simulations start with the no-blowing ca se s , illustrating the flow charac- 
teristics in pre-stall and post-stall conditions. Tangential blowing is implemented 
under the same free stream conditions, and the effects of blowing on the vortical 
flow Add are discussed. 

4.2 Multiple-Grid Concept and Grid Generation 

A proper arrangement of grid distribution contributes to a successful numerical 
simulation. Far solving the tangential leading edge blowing problem, a desirable 
surface-grid arrangement on a cross-flow plane is shown in figure 41, where the 
grid paints are clustered near and at the jet slot. 

The use of a single-zone, as described in Chapter 3, poses a major difficulty : the 
grid points tend to cross one anothe r near the jet-exit due to both the large number 


Once the grid distribution on this inner boundary is specified. Grid 1 is generated 
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through a two dimensional hyperbolic grid-generation method [35], aa described in 
Chapter 3. 

The second zone (Grid 2) fills the upper half a t the open which es- 

sentially co ver s the jet regno and co-flowing region (where the jet and outer flow 
merge), as illustrated in figure 4.3. Similarly, Grid 3 is constructed by filling the 
lower half of the open domain such that Grid 2 and Grid 3 wrap entirely around 
the delta wing. Both Grid 2 and Grid 3 are generated algebraically under the con- 
straint of the orthogonality condition normal to the walL An exponential stretching 
function is used in these two zones to duster the grid points near the walL At each 
cross-flow plane, a common boundary line at the leading edge is shared between 
Grid 2 and Grid 3. The number of grid paints of each zone that lie on a com- 
mon boundary line is usually different, and their locations do not coincide with one 
another. 


The final step is to redistribute the grid points along each radial curve in zone 1 
such that Grid 1 starts from the last interior points of Grids 2 and 3. In other 
words, there are 2 grid points in Grid 1 that coincide with those in Grids 2 and 3. 
The spacing of this overlap region (figure 4ul) is determined as follows: 


«i 


+ <s 
2 


(4.1) 


where the a v er a ged grid spacing of Grid 2 (e a ) and Grid 3 (e 3 ) just outside of the 
overlap region is used. This is to ensure a smooth transition of grid spacing between 
zones. An exponential stretching function is also used in Grid 1 to extend the grid 
points to the outer boundary, with minimum spacing based on the length of the 
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overlap region («t). 

The final grid arrangement ia idmilar to that shown in Chapter 3 (figure 3.3) 
except that the number of zones and the total number of grid paints are different. 
The breakdown of the grid paints for each zone is summarised as follows: Grid 1, 
35x35x33 - 40,425; Grid 2, 35x18x31 * 19, 530 and Grid 3, 35x35x25 * 15,750, 
which results in a total number of 75,705 grid points. For each Grid, the numbers 
of grid points correspond to the {, *) and ( directions respectively. 


4.3 Boundary and Jet-Exit Conditions 

The boundary conditions used in this computation are essentially the same as those 
described in Chapter 3 with the exception of interface boundary conditions, en- 
countered in the multiple-grid scheme. Since grid points coincide with one another 
in the overlap region as described in the previous section, the interface boundary 
conditions for Grid 1 are obtained directly from those of the interior points of Grid 
2 and Grid 3. Conversely, the interface boundary conditions for Grid 2 and Grid 3 
are determined from Grid 1. A linear interpolation is used at the boundary 

between Grids 2 and 3. Since the number of grid points and their locations in 
these two zones are not the same at the interface boundary, the conservative flow 
variables ($) are fitted with a cubic spline before the interpolation is performed. 
The boundary conditions are summarized as follows: 

(1) Free stream values are specified at the upstream and circ u mferential bound- 
aries and remain fixed. 
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(2) Extrapolated values are used at the downstream boundary except for e, which 
is calculated based on the free stream static pressure for the present subsonic 
calculations. 

(3) Symmetry conditions are imposed on the symmetry plane. 

(4) No-slip, adiabatic wall and zero pressure gradient conditions are imposed at 
the wall, as discussed in section 3.3, except for the jet region when Mowing is 
implemented. 

It should be noted that the presence of the wall is felt by Grid 1 only through the 
interface overlap regions with Grid 2 and Grid 3. 

The blowing momentum coefficient for the conical m omentum distribution, for 
a jet exiting with the free-stream density is defined as follows 



where is the jet-exit velocity, A is the local jet-slot width and s is the local 
span. Since conical flow conditions were intended in the experiments [18, 19], both 
the external dimension of the delta wing and the jet-slot width increase linearly 
from the apex, with (£) having a constant value of 0.0046. The jet-exit boundary 
conditions are summarized as follows: 

(1) The local static pressure is used at the jet-exit. This is approximated by 

the static p re ss ure value at the nearest point in the span wise direction 
from the lending edge. 

(2) The free stream value of density is specified at the jet exit. 
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(3) The jet-exit velocity distribution is prescribed as a power function given by 


V,WV M i 


-(¥)T 


(4.3) 


which is modified from the ex pr es si o n for an ellipse (0 = \). A distribution 


was chosen with 0 * 0.1 to simulate a turbulent wall jet velocity profile, as 
shown in figure 4.4. The maxiimim velocity magnitude is determined accord- 
ing to a specified from Eq.(4.2). 


4.4 Results and Physical Interpretation 

4.4.1 No Blowing 

A 60*-sweep, rounded leading edge delta wing is used in the computations. The 
delta wing is modeled to have conical geometry up to 80 % of Co with a maximum 
thickness of 12.5 % of at z/co = 0.8. The aft part of the wing consists of an 
afterbody whose thickness varies linearly toward the trailing edge, as shown in figure 
4.5, to simulate a sharp trailing edge wing. This wing has the same proportions 
as the one used in Wood's experiments [18, 19] except that the experimental model 
didn't have the afterbody, thus ending in a blunt trailing edge, as shown in figure 
1.4. 

The experimentally prescribed values [18, 19] of a » 30* and a =» 40*. corre- 
sponding to pre-stall and post-stall conditions respectively, are used in the present 
no-blowing simulation, so that the Sow characteristics under these conditions can 
be distinguished. 
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Figures 46(a), (b) show "particle traces” from particles released near the leading 
edge for 30* angle of attack, in plan and side views. The formation of a well 
organized vortical flow Add is observed in these figures. A conical jet-slot dimension 
is simulated where the slot width increases linearly from about 20% to 80% (where 
the maximum body thickness occurs) of the root chord. It is interesting to note that 
some of the particles just in front of the afterbody are trapped in a redrculatory 
motion near the jet-slot, where solid wall boundary conditions are implemented. 
Releasing the particles from the same locations at a * 40*, as shown in figures 
4.7(a), (b), a chaotic behavior of the flow is revealed above the wing, indicating a 
reverse flow region associated with vortex breakdown. This representation of vortex 
breakdown can often be found in recent pu b lications [31, 54]. From observation of 
the cross-flow velocity profiles at a at 40*, (to be shown in the next section), the 
flow separates near the leading edge just ahead of the jet-slot, so that particles are 
no longer trapped. 

Since the "particle traces” depend on the particle release location, these particle 
paths can only qualitatively represent the region where vortex breakdown occurs. 
A more precise repr e sen tation is to plot the cross-sectional reverse-flow velocity 
( negative- u velocity) contours at a number of selected stream wise locations, as show 
in figure 48. These streamwise stations are chosen arbitrarily at x/co w .17, .29, .42, 
.54, .67, .79, and .92. For a = 30*, the r eve r se flow bubble does not appear (not 
plotted), whereas at a = 40* large regions of reverse flow are observed indicating 
extensive vortex breakdown. In addition, the parameter u at each cross-flow plane 
is defined as the ratio of the maximum reverse-flow (negative-u) velocity to the free 
stream velocity 
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S- aegsfe (4.4) 

On each crow-flow plane, the outer contour haa aero negative-u velocity, which can 
be interpreted aa the reverse flow bubble boundary. The values of u indicated in 
figure 4.8 refer to the innermoat contour of the correaponding station. 

Another phenomenon associated with vortex breakdown is the rapid expansion 
of the vortex core region. This can be demonstrated by means of total pressure 
contours, aa shown in figure 49, for the station at x/e® * 0.36. For a » 30*, the 
vortex is highly concentrated with a strong low in total pressure in the vortex core 
region. For a = 40*, the vortex appears diffused, the low in total pressure is lew 
significant and extends over a much larger region. 

Figure 410 (a) shows the comparison of computed span wise premure distribu- 
tion for a ss 30*, x/e® = .36 with that from experiments [18, 19]. Agreement is 
good in the sense that the predicted and measured location and strength of the 
vortex-induced suction peak are similar. The discrepancy near the root chord may 
be attributed to the displacement thicknem associated with the formation of the 
boundary layer on the wind tunnel side wall, which may have translated in an out- 
board displacement of the pressure curve. Another factor could be the lack of grid 
resolution in that area. 

The comparison for a = 40* at the same chord location is shown in figure 
410 (b), where good agreement with experiments is also obtained. In this case, 
a nearly fiat pressure distribution on the wing upper surface signifies a weak and 
diffused vortical flow, another indication that vortex breakdown has occurred. 
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Figures 4.11 (a) and (b) show the computed near-surface streamline patterns 
for a as 30* and a am 40". At a at 30", a typical vortical flow pattern on the upper 
surface similar to that described in Chapter 3 is observed. At a » 40", however, 
the advene pressure gradient in the longitudinal direction associated with vortex 
breakdown, accompanied by the diminiahing influence of the vortical flow, causes 
reverse flow on the wing surface. The afterbody separation due to the rapid change 
in the slope of the wing surface at 80% of the root chord is apparent in both cases. 

The computed normal force coefficient versus angle of attack for no blowing is 
shown in figure 4.12, where C w_ T r occurs approximately at a * 30*. The decrease 
in normal force at a — 40", due to massive vortex breakdown is evident. 


4.4.2 Tangential Blowing 

The same delta wing and grid arrangement are used for the tangential blowing 
simulations. The conical jet slot geometry extends from 17% to 79% of the root 
chord with a constant jet-slot width (h) to local wing span (s) ratio, h = 0.0047. 
The same angles of attack (a = 30* and 40"), are used in the blowing calculations 
as in the no-blowing cases, so that the effects of tangential blowing can be identified 
through parallel comparison. 

At a at 30" with C„ a 0.05, figure 4.13(a) shows the side and plan views of 
"particle traces” emanating from the same locations as those in the no-blowing case 
(figure 4.6(a)). Very little difference is observed between the blowing and no-blowing 
cases, except that the particle traces appear to reduce in extent above the wing due 
to blowing and particles are not trapped near the jet-slot as in the no-blowing case. 



Figure 4.17 shows the comparison of the computed suz&ce pressure distribution 
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with Wood's experimental results at x/co * 0.36 for a * 30* and 40*. The agree- 
ment is good for both angle-of-attack cases. However, the predicted suction peak 
associated with the jet is not as high as that obtained experimentally. This dis- 
crepancy may be attributed to the lack of grid resolution near the jet region. Grid 
refinement studies were not conducted due to limitations of computer memory. 

Figure 4.18 shows the computed span wise surface pressure distribution at x/co = 
0.36 for a s 30*, with and without blowing, to illustrate the effects of blowing in 
pre-stall condition. It is observed that blowing tends to reduce the vortex-induced 
pressure peak while increasing the suction peak around the leading edge. This 
indicates that the delayed separation of the cross-flow boundary layer modifies the 
strength and position of the resultant vortex as the flow Arid adjusts itself to a 
new equilibrium state. Figure 4.19 illustrates the changes of the primary separation 
point in the cross-flow plane, with and without blowing. 

Figure 420 illustrates the effect of blowing momentum on the spanwise pressure 
distribution at x/cq = 0.36 for a a 40*. The effect of blowing is quite apparent: 
blowing re-establishes the vortical flow, as indicated by the vortex induced suc- 
tion peak, and creates an additional suction peak around the leading edge. These 
phenomena may be attributed to the reduction of vortex strength brought about 
by altering the cross-flow separation point through the wall jet, as shown in figure 
421. In the absence of Mowing, the vortex would burst. As a result of the re- 
establishment of orderly vortical flow, the integrated value of the surface pressure 
increases, which translates into higher normal force an the wing. 


The effect of blowing an the computed near-surface flow patterns on the upper 
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surface of the delta wing for a * 30* and 40* are shown in figures 4.22 (a) and 
(b). At a > 30*, the surface flow pattern is to that in the no-blowing ease 

(figure 4.11(a)), except that the secondary separation line appears to move inboard 
slightly indicating an inward motion of the vortex due to blowing. For a = 40*, 
the recovery of the vortical flow pattern due to blowing is quite obvious as 

compared to the noblowing case (figure 4.11(b)). The re verse flow region due to 
vortex breakdown in the noblowing case has been removed, an indication of the 
restoration of orderly vortical flow. Due to the rapid change in body geometry and 
the finite extent of the leading edge jet (17% - 79% of Co), the afterbody separation 
still remains, as shown in figures 4.22 (a) and (b). 

Figure 4.23 illustrates the computed normal force coefficient as a function of 
the angle of attack for noblowing and blowing with a constant blowing momentum 
coefficient of 0.05. Three observations can be made about these results : 

(1) Tangential blowing delays vortex breakdown at high angles of attack, thereby 

increasing the stall angle or the angle of attack at which Cm occurs. 

(2) Tangential blowing restores the vortical flow in the post-stall conditions, in- 
creasing the maximum normal force an the wing. 

(3) In pre-stall conditions, tangential blowing tends to reduce the pressure peak 
associated with the vortex and to increase the suction peak around the leading 
edge such that the integrated value of the surface pr es sure remains about the 
same. In other words, these two pressure peaks ofiset one another resulting 
in a nearly constant pressure-induced lift. 
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(a) Side view. 



(4) Plan view. 


Figure 4.6: Illustration of vortex structure with particle traces Cor a = 30*, C„ = 0. 
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(a) Side view. 



(b) Plan view. 


Figure 414: Illustration of vortex structure with particle traces, for a s 40*. 
C„ = 0.038. 
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Figure 4.15: Reverse-flow (negative-u) velocity contours at selected chord wise loca- 
tions, x/cq = .17, .29, .42, .54, .67, .79, .92, for a = 40*, C 0 = 0.038. 


(6) a = 40*, C„ = 0.038. 


Figure 4.16: Normalized total pressure contours at x/cq 
C M = 0.05 and a = 40", C M = 0.038. 
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Figure 418: Computed span wise surface pressure distributions for a = 30*, 
x/cq = 0.36, with C M = 0 and C„ = 0.05. 
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Chapter 5 


Conclusions and 
Recommendations 


5.1 Conclusions 

Numerical simulations of spanwise and tangential leading edge blowing on delta 
wings have been successfully carried out by using Computational Fluid Dynamics 
techniques for solving the three-dimensional Thin-Layer Navier-Stokes equations. 
This is a substantial departure from the previous theoretical studies, where the 
assumption of a conical flow field was a requirement for solving the blowing problem. 
A single grid is used for the spanwise blowing cases, whereas a multiple-grid zonal 
approach is utilized in the tangential blowing calculations to acc om odate the jet-slot 
geometry. In both cases, the leading edge jet is simulated by defining a permeable 
boundary, corresponding to the jet slot, where suitable boundary conditions are 
implemented. The algebraic turbulence model for a wall jet is inc o rporated to 
accurately predict the location of the primary separation point in the tangential 
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blowing calculations. For both blowing schemes, the numerical results axe in good 
general a g re emen t with data from wind tunnel tests. An understanding of the 
vortical flow field structure over a wide range of angle of attack, with and without 
blowing, was facilitated by detailed graphical repr e s entations. 

It is found that spanwise blowing is useful at low to moderate angles of attack, 
where an increase in lift is obtainable through the enhancement of the strength of 
the leading edge vortices. The main conclusions for spanwise leading edge blowing 
are summarized as follows : 

(1) Spanwise blowing is mainly an inviscid phenomenon, as suggested by previous 
in viscid calculations and verified by the present thin-layer approximations. 

(2) Blowing increases the lift by magnifying the size and strength of the leading 
edge vortices at a given angle of attack. As a consequence, the increase in lift 
is primarily due to the non-linear vortical contribution. 

(3) Blowing displaces the flow outboard corresponding to an effective enlargement 
of the wing-span. 

(4) Blowing causes vortex breakdown at a lower angle of attack. 

The use of tangential leading edge blowing is found to be useful at high angles 
of attack for delaying vortex breakdown and restoring orderly vortical flow. The 
effects of tangential blowing are summarized as follows : 

(1) Tangential blowing is mainly a viscous phen om enon since the tangential jet 
delays primary separation. 
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(2) If the vortical flow is orderly (i.e. pre-stall regime), blowing tends to reduce 
the vortex induced p res s ure peak and to increase the jet induced suction peak 
around the leading edge such that the integrated value ai surface pressure 
remains about the same. 

(3) Blowing delays vortex breakdown. 

(4) Blowing preserves orderly vortical flow to higher angle of attack, leading to 

an increase of Cm and stall angle. 


In summary, spanwise and tangential leading edge blowing can be used, within 
different ranges of angle of attack, to control the vortical flow over delta wings or 
delta-like lifting bodies. The physical phenomena for these two kinds of blowing are 
quite different in the sense that the effects of spanwise blowing are dominated by 
inertia forces while the changes in the structure of the flow field due to tangential 
blowing are brought about by viscous effects. It appears that spanwise blowing is 
useful at low to moderate angles of attack to gain lift. On the other hand, tangential 
blowing is useful at high angles of attack to delay vortex breakdown and to restore 
the orderly vortical flow field. The use of either one of the blowing schemes can be 
beneficial depending on the operational envelope and angle of attack regime of the 
aircraft. 


5.2 Recommendations 


The first recommendation is to do grid refinement studies on both blowing schemes. 
The lack of sufficient grid resolution may, in part, provides the discrepancy between 
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computational result* and experimental measurements. One approach is to utilize 
the Solid State Device (SSD) on the Cray X-MP for the multiple-grid zonal program 
such that zonal information not currently in use could be stored in the SSD. In this 
way, less demand for on-core memory would allow an increase in the number of 
grid paints. An easier approach would be to use the Cray-2 machine, which has 
substantially more memory than the Cray X-MP. 

The second recommendation is to simulate the blowing schemes by solving the 
full Navier-Stokes equations. For high angles of attack, massive separation extends 
over a large region where the coarse grid spacing may not be sufficient to resolve 
high flow gradients associated with a wake structure. Inaccuracies stemming horn 
the thin layer approximations may be another cause for discrepancy between com- 
putation and experiment. 

Concerning future work, two additional ffwnnwn^tinni are made. 

(1) Asymmetric blowing conditions should be considered in order to study blowing 
as a means of roll control 

(2) The numerical results should be used to formulate simple scaling laws relat- 
ing the pre-stall and post-stall aerodynamics to angle cl attack, aspect ratio 
and blowing momentum coefficient. The correlation of these parameters from 
numerical results would allow for a substantial reduction of effort in an ex- 
perimental program. 




Appendix A 
Jacobian Matrices 

The three-dimensional, inviscid flux Jacobian matrices in generalized curvilinear 
coordinates are given by ( [37]) 

K t K, 

Ktf* — u8 K t + 8 — K,(7 - 2) U 

A, B OT C — Kjj? — t l9 KgV — Kf (7 — 1) U 

K,<t 7 — rv9 K X W — /C, (7 — l)u 

9 (2<P - ^ep- 1 ) k, (yep- 1 - 4> 2 ) - (7 - 1 ) uO 

#Cj 0 

«,«!-*, ( 7 - 1)0 K,u-K M (y-l)w *,( 7-1) 

« f + tf-a f (7-2)t; (7 - l)w ( A>1 ) 

(7 -l)o K t + 9 - K M ( 7-2)10 *» (7 — 1) 

** (7e/» -1 - ^ 2 ) - (7 - 1) *« (l e P~ l - 4?) - {y - l)w8 K t + y9 

where 

9 SB «,U + + K,U7 
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4 


4 >* - ;(?-!)(«* + »’ + «’) 


with K»{,T|or({orA,lott, respectively. 

The viscous flux Jacobian matrix [36] is 

0 0 0 0 0 

mn <*\6<(p- x ) a A (/»■*) «s*< (P~') 0 

iCr*»lmst ai6((p~ l ) a<6 ( (p~ l ) a, S < (p~ l ) 0 

rrtu as6 ( (p- 1 ) a,S((p- 1 ) as S ( (p~ l ) 0 


where 


m*i 


mu 


ntu 


mu a 0 «<(/>"*) 


m n * (u/p) - ajS ( (v/p) - a^ c (w/p) 

mu * -a 3 S ( (u/p) - a 4 $ ( (w/p) - a t 6 c {w/p) 
m« » — aj^< (u/p) - a,^ (o/p) - as$c (»/p) 
m M * -a l S ( {u i /p)-a r S c (2uv/p)-a 3 S((2uto/p) 
-a 4 S ( (o’/p) - a*# c (w’/p) - a 5 * c (2w>/p) 
-<*<£( (e/p’) - ao^c [(u* + o’ + to*) /p] 
m*a * -m n - a«^ c (u/p) 
mu = -mn -aoS ( (v/p) 
m M = —nisi — ao^( (w/p) 

as - 7«i*r-‘(£ + $ + C?), «4-#*(cJ + 5< t a + C. a ) 

<*\ = p + cj + c?) » <** = jC»c* 

«a = §«,. «s-M(c2 + < + |c. a ) 

a, = 


(A.2) 





Appendix B 


Similarity Transformation 
Matrices 


The similarity transformation matrix [37] for the Jacobian matrices A, B or C is 


T, = 


** 

k a u 

** v + k,p 
k,w-k,p 


«*tt - KmP 
k,v 

kyW + k x p 


[ + />(**» - *,»)] [*, p£iy + P( k * v ~ *«»)] 


a,u + kfp 

k,v - k?p 
k,w 


a 

a (u + kga) 
a(e + k^a) 
a (u> + k t a) 


a 

a (u + kgo) 
a(v + k^a) 
a(w + k,a) 


[«.i£y + />(*,» -*.»)] « [?£§ + «*] a [l£S + 0 *]J 


(B.1) 
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where 

*« 
a 

** 

9 

K 

The analytical expression for the inverse of T, i* 

[*» (l - £) - P~ l ( k,v - *,«;)] *, (7 - 1 ) uaT 2 

[*» (l “ $r) “ P' 1 (*»w ~ *.«)] [~k,p~ l + S( 7 ~ 1) «a- J ] 

[*« (l - S) - P~ l (*»« “ *«»)] [-Kjp -1 + *,(7 - 1) «<*"*) 

0 [kga — (7 — 1) u ] 

Pfat + ai) —0 [k,a + (7 + 1) u] 

[k,pT x + k, (7 - 1 ) t>a“ J ] + *« (7 - 1) wo' 1 ] -*« (7 - 1) <*"* 
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[-k*P~ l + *» (7 “ 1) *» (7 - 1) w«T 2 -k x (7 - 1) a"* 

0(**a-(7~ 1)»] (7- 1)»] 0(7-1) 

~"P [*»« + (7 ~ 1 ) ®1 —0 [*» a + (7 ~ 1 ) H 0 ( 7 - 1 ) 

(B-2) 





where 0 = 1/ y/2pa 
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The analytical expressions for {Tj~ l T n ), (T^ l T ( ) used in Eq. (2.39) and their inverse 
matrices (T^ l T ( ), (I^T,) can be written in a simple form. 


37 'T, 


where 


rrn 

mj 

m 3 
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-mj 

n»i 

m 4 

#*m 3 

-ltm s 
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m, 
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Calculation of the Jet Momentum 


The jet m om e ntum (J) in the absence of compressibility can be exp r e ss ed as follows 


J = p J J V if dydx 


(C.1) 


The expression for the velocity distribution for a mull trapezoidal cell, as shown 
in figure C.l, can be written as 


v ^-W> v+v ' 


(C-2) 


where Vi, V 2 , Vs and V« are the velocity magnitudes at the cell corners, with 


n(*) 

Vs(*) 

Ay(*) 


Ax 

Va-V* 

Ax 

U-lx 

Ax 


* + V, 

* + V, 

x + Li 


(C.3) 
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APPENDIX C. CALCULATION OF THE JET MOMENTUM 


Integrating Eq. (C.2) at a constant * yields 


- ^^[v.’ + v,v.+v>] 


Substituting the e x pr es si ons lor Vi(x), V*(x) and Ay(r) from (A.3) into (A.4) and 
integrating along the x-axis yields 


r&M ft±*(a) 

J « PI J V*(x,y) dydx 
jo jo 


£jj>F^ + (PJ5? + P0)^ + (QE + FJl)^ + JtEAx (C.5) 


where 
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